In-plane magnetic anisotropy of Fe atoms on Bi 2 Se3(lll) 
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The robustness of the gapless topological surface state hosted by a 3D topological insulator against 
perturbations of magnetic origin has been the focus of recent investigations. We present a compre- 
hensive study of the magnetic properties of Fe impurities on a prototypical 3D topological insulator 
Bi 2 Se3 using local low temperature scanning tunneling microscopy and integral x-ray magnetic cir- 
cular dichroism techniques. Single Fe adatoms on the Bi2Se3 surface, in the coverage range ~ 1% are 
heavily relaxed into the surface and exhibit a magnetic easy axis within the surface-plane, contrary 
to what was assumed in recent investigations on the opening of a gap. Using ab initio approaches, we 
demonstrate that an in-plane easy axis arises from the combination of the crystal field and dynamic 
hybridization effects. 



Topological insulators (TI) have demanded heavy in- 
terest from the scientific community as a new class of 
materials illuminating fascinating yet exotic physics and 
offering large potential for applications in the field of 
spintronics [1 . TI host a gapless topological surface state 
(TSS) which exhibits a Dirac-cone like dispersion. How- 
ever, unlike in the case of graphene, the Dirac cone is 
located in the center of the Brillouin-zone and the spin 
and momentum degrees of freedom are locked. The lat- 
ter so-called topological quality of the TSS is protected 
by time-reversal symmetry which leads to a variety of 
interesting effects. For example, these materials may be 
a forum for Majorana fermions [2], a topological mag- 
netoelectric effect [3 , and a quantized anomalous Hall 
effect 0]. 

The locking of both spin and momentum degrees of 
freedom leads to the suppression of 180° elastic back- 
scattering of TSS electrons in the absence of spin- flip pro- 
cesses. The robustness of these "topologically protected" 
processes, when introducing impurities which break time 
reversal symmetry, is of critical importance for spin- 
based transport in such materials. It has been suggested 
that the interaction between magnetic impurities and the 
topological state can cause an opening of an energy gap at 
the Dirac point (DP), provided that the magnetic order 
is oriented normal to the surface plane [5H7]. Neverthe- 
less, the stability of the TSS against local magnetic per- 
turbations is currently under heavy debate [8HT0], since 
little is known about the fundamental interface effects of 
magnetic entities with TI surfaces as well as the static 
magnetic properties of impurities in these materials. 

In this work we present a comprehensive study of the 
magnetic properties of iron impurities on a prototypi- 



cal TI Bi2Se3 using local scanning tunneling microscopy 
(STM) and integral x-ray magnetic circular dichroism 
(XMCD) techniques. We show that the multiplet struc- 
ture visible in x-ray absorption spectra reflect a high-spin 
state of the adsorbed Fe impurities, which are confirmed 
to bind at hollow sites of the surface lattice, under the 
influence of trigonal crystal fields given in hollow site 
positions of the Bi2Se3 surface. In contrary to recent re- 
ports a significant in-plane magnetic easy axis direction 
is observed in the XMCD. This observation is supported 
by our DFT calculations which reveal the possibility of a 
preferred in-plane easy axis when considering both crys- 
tal field and dynamic hybridization effects. 

Both local scanning tunneling spectroscopy (STS) and 
integral XMCD measurements were carried out with one 
and the same Bi2Se3 crystal used in Ref. [11]. Clean 
Bi 2 Se3(lll) surfaces were obtained by in situ cleaving at 
room temperature under ultra-high vacuum conditions 
followed by a fast cool-down to temperatures of about 
T = 10 K within 15 minutes. In order to prevent sur- 
face diffusion, Fe was deposited at lowest temperatures 
onto Bi2Se3(lll) using an e-beam evaporator. For the 
case of x-ray measurements the general procedure of in 
situ metal atom deposition and coverage calibration is 
described in detail in Ref. [T2] . 

STS was performed at a temperature of T — 0.3 K 
using a home-built ultra-high vacuum STM with a tung- 
sten tip [13 . Figure [TJ a) shows « 1% of a monolayer 
(ML) Fe deposited on a clean Bi2Se3 surface imaged with 
constant-current mode. Single Fe atoms are located with 
an almost equal ratio in two different hollow site posi- 
tions with respect to the underlying hexagonal lattice 
of Se surface atoms. High resolution images (Fig. [lib)) 
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FIG. 1: (a) STM constant-current image of 1%ML of Fe 
atoms adsorbed onto a cleaved Bi2Se3 surface; 7 s tab = 500 pA, 
V^stab = +100 mV. (b) High resolution image of two Fe 
atoms at different binding sites (A,B); i s tab = 600 pA, 
Ktab = —100 mV. (c) Tunneling spectra taken on two Fe 
atoms at the indicated binding sites and on a substrate spot 
of the Fe covered surface. The gray bar refers to the DP 
of the uncovered Bi2Se3 surface as confirmed by probing the 
zeroth Landau level in high magnetic fields. Inset: same as 
(c) after subtracting the substrate spectrum; I sta b = 600 pA, 
Vitab = -100 mV, V mod = 3 mV; / mod = 4.1 kHz. 



show that the atoms on the two binding sites, labeled A 
and B, clearly differ in apparent height and shape. The 
apparent height of 0.5 A and 0.3 A for A and B, respec- 
tively, is small compared to e.g. Fe on Cu(lll) (~ 1 A, 
[14]). A closer inspection reveals that the apparent shape 
of the Fe and of the halo (charge density) surrounding the 
atom depends on the binding site. Type A atoms show 
a circular protrusion around the center of the Fe atom 
while type B atoms show a triangular protrusion. The 
surrounded halo pattern is quite complex but unique for 
each binding site. Infrequent hopping events of a single 
Fe atom between A and B stacking positions indicate a 
transformation of atom shape and halo between the two 
distinguishable atom types and binding sites. This proves 
that these properties are not due to underlying surface 
defects but are purely electronic effects resulting from 
the bonding between the Fe atom and the atomically flat 
Bi2Se3 surface. 

Tunneling spectra (dI/dV(V)) that were taken in dif- 
ferent magnetic fields perpendicular to the surface on the 
bare Bi2Se3 surface prior to Fe deposition (not shown) ex- 
hibit Landau levels (LLs) with the expected energy scal- 
ing confirming the Dirac cone dispersion of the topologi- 
cal surface state [15j [16] as well as the position of the DP 
which is measured in a window around £?dp = -0.31 eV. 
This is consistent with ARPES measurements of the band 
structure of the same crystal [11 . This DP binding en- 
ergy is typical for naturally grown Bi2Se3 crystals, which 
are known to be electron doped due to the presence of Se 
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FIG. 2: X-ray study of Fe impurities (coverage 1%ML) on 
Bi 2 Se 3 at T = 10 K: (a) Magnetization curves M(B) for di- 
rections = 0°, 70° in red and blue, respectively. The dotted 
lines represent thermodynamic fits with a moment m — \\iB 
and K — 1.6 meV (see text), (b) XMCD signals in a field 
of B — 5 T normalized to the XMCD peak intensity at the 
L2-edge. (c) and (d): XAS (red, upper panel) and XMCD 
(black, lower panel) measured in a field of B = 5 T for di- 
rections — 0° (c) and 70° (d). Green spectra refer to CF 
calculations with D a = 0.19 eV and D T = 0.05 eV. 



vacancies observable in STM images (Fig. Qa)) [TTt fT8] . 
Upon deposition of Fe (Fig.JIJc)), the substrate spectral 
intensity reveals a minimum that is shifted downwards 
with respect to the uncovered surface by about 50 meV. 
Since, for the uncovered surface, this minimum coincides 
with the DP, we conclude that Fe probably acts as a 
donor and the DP is accordingly shifted downwards as 
seen by recent ARPES measurements [9 j. Most notably, 
a resonance can be seen in spectra taken on the Fe atoms. 
It appears at voltages slightly above the DP and has dif- 
fering intensities, peak positions, and linewidths depend- 
ing on the binding site. Substrate corrected spectra (in- 
set) prove that the resonance has a considerably stronger 
intensity and is shifted closer to the DP for the type B 
Fe atoms. Similar resonances have been observed on the 
Se vacancies [T7J [18] and were attributed to Coulomb 
scattering of the topological surface state by the impu- 
rity H [18]. 

For a magnetic characterization, x-ray absorption spec- 
tra (XAS) at the Fe Ls^ absorption edges (between 
690 eV and 740 eV) were recorded in the surface sen- 
sitive total electron yield mode for different angles be- 
tween the Bi2Se3 surface normal and the x-ray beam di- 
rection. Magnetic fields B were applied collinear to the 
beam as shown in the inset of Fig.|2|a). The experimen- 
tal timescales for measuring spectra are on the order of 



100 s. Before and after the magnetic characterization, 
the samples were carefully checked for oxygen contami- 
nations using the XAS at the 0-1 s absorption edge. 

Typical non-dichroic XAS defined as the average of 
XAS for right (<r + ) and left (<J~) circular polarization are 
shown after Bi2Se3 background subtraction in the upper 
panel of Fig. |2|c) and (d). When B is applied, a pro- 
nounced XMCD signal defined as (a + — a~) is observed, 
which reflects the spin and orbital polarization in the 
Fe 3d states. As a first approximation the XMCD signal 
is proportional to the average Fe magnetization M when 
it is normalized to the L3 peak intensity in the XAS signal 
as shown in the lower panels of Fig. [2^c) and (d) for fields 
of B = 5 T. More precisely, XMCD/XAS|l 3 oc E • M, i.e. 
the XMCD signal normalized to the total XAS intensity 
is proportional to the projection of the averaged Fe mo- 
ments < m > on the photon propagation direction E. 
Comparing the lower panels of Fig. [2^c) and (d) we thus 
conclude that Fe atoms prefer to be magnetized in the 
= 70° direction as compared to the out-of-plane case 
where (9 = 0°. 

This is supported by the significant anisotropy in the 
orbital moment tul visible in Fig.[2^b): the XMCD spec- 
tra are normalized to their peak intensity at the L2 edge 
(721.0 eV) and smaller intensities in the polar geome- 
try at the L3 edge (708.0 eV) indicate considerably re- 
duced orbital moments tul with respect to the direction 
= 70°. By measuring the normalized XMCD intensity 
at the L 3 edge for different magnetic fields it is possi- 
ble to obtain magnetization curves M(B). The curves 
M(B) in Fig. [2|a) at = 0°,70° show smaller values 
along = 0° again proving a preferential orientation of 
the magnetization in the (111) plane of Bi2Se3. 

In the following we want to discuss the experimen- 
tally observed magnetic anisotropy of the Fe ad atoms. 
Let us first look at the XAS lineshape in Fig. [2jc) and 
(d), which reflects the electronic configuration and pos- 
sible crystal field (CF) effects from the Bi2Se3 surface. 
A pronounced double-peak structure is observed at the 
L2 edge as well as a shoulder at the high-energy side 
of the L 3 edge (see arrows). CF multiplet calculations, 
including electron repulsion and spin-orbit coupling [19] 
can reproduce the main features assuming a purely trig- 
onal CF of D3^ symmetry: the green curves represent 
spectra derived for Fe d 6 impurities in a CF defined 
by the parameters D a = 0.19 eV and D T = 0.05 eV. 
The double-peak structure in the XAS of the L2 edge at 
= 0° as well as the shoulder at the L3 edge for both 
= 0° and 70° are reproduced. The measured spectra 
are broadened with respect to the model, which points 
towards a significant coupling of d-states to delocalized 
bands. From our CF calculations also considering Boltz- 
mann statistics at T = 10 K, we derive anisotropic or- 
bital moments with Am^ =0.13 \±b in the two measur- 
ing geometries supporting the experimentally observed 
quenching of the orbital moments in the = 0° direc- 
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FIG. 3: (a), (b): Adsorption geometries of Fe (red) on Bi2Se3 
(Bi: gray, Se: yellow). Fe can adsorb at an fee (a) or an hep 
hollow site (b) . The real and imaginary parts of the hybridiza- 
tion functions (ReAu(E) + e% and ImAn(E)) are given in 
(c,d) and (e,f), respectively, and quantify the interaction of 
the Fe 3d-electrons with their environment. Au(E) is shown 
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tion under the influence of the in-plane CFs. Assuming 
a moment of \m\ = 4 \±b we can fit the experimental 
magnetization curves in Fig. [2|a) with a thermodynamic 
model which includes the Zeeman term and an anisotropy 
term Ea = K - (z - m / \m\) 2 where K is the magnetic 
anisotropy per Fe atom and z the normal of the (111) 
plane. The fit leads to a significant in-plane magnetic 
anisotropy of K = +1.6 meV per atom. 

To understand the physical mechanisms behind the ex- 
perimentally observed magnetic properties of the system 
we resort to ab initio calculations and derive a quantum 
impurity model describing the Fe adatoms on the Bi 2 Se3 
surface. We performed density functional (DFT) calcula- 
tions of Fe adatoms on the Bi2Se3 surface which yield the 
relaxed adsorption geometries, the magnetic moments, 
and the hybridization functions [20l [21] [27] . 

Our calculations reveal two possible adsorption sites 
for an Fe adatom at an fee hollow (f-) or an hep hollow 
(h-) site [Fig. [3] (a,b)]. While both adsorption geome- 
tries are stable, we find that the f-site adsorption yields 
0.07 eV lower total energy than the h-site. Independent 
of the adsorption site, GGA predicts a Fe 3<i-occupancy 
of Nd = 6.3. At both sites, we find a strong relaxation of 
the adatom into the surface with equilibrium positions of 
the Fe adatom virtually at the same height (±0.1 A) as 
the surrounding Se atoms. This is well in line with the 
STM topography (Fig. |TJl) and the observed broadening 
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in the x-ray absorption spectra. It has important conse- 
quences for the electronic and magnetic properties of the 
Fe adatoms. 

The energy dependent coupling strength of the impu- 
rity 3<i-orbitals and the surrounding crystal is quantified 
by the hybridization functions Aij(E) where z, j = (m, a) 
denote combined orbital and spin quantum numbers of 
the Fe 3d electrons and E is the energy. These hybridiza- 
tion functions can be interpreted as energy dependent 
complex valued potentials acting on the Fe 3d orbitals 
and fully characterize their interaction with the substrate 
as regards local observables (see Ref. [22] and the sup- 
plement for further details). For Fe in the h-site posi- 
tion (Fig. [3^f)), we find the most pronounced peaks in 
ImAu(E) for the in-plane d-orbitals (d xyi d x 2_ y 2 span- 
ning the subspace of orbital angular momentum quan- 
tum number \l z \ = 2) at E « — 2 eV and for the out- 
of-plane oriented d z 2 orbital above the Fermi level. For 
Fe at the energetically most favorable f-site (Fig. [3^e)), 
ImA^(£) is clearly dominated by peaks at E w — 2 eV 
in the d xy and d x 2_ y 2 orbitals. The hybridization of the 
d z 2 orbital is smaller than at the h-site, because there is 
no Bi atom directly beneath the Fe atom. Consequently, 
at both adsorption sites, the in-plane oriented d xy and 
d x 2 _ y 2 -orbitals hybridize strongly with electronic states 
from the surrounding crystal. The relaxation of the Fe 
adatoms into the Bi 2 Se3 surface facilitates the strength 
of this particular hybridization channel. 

Re An (E) can be interpreted as the energy dependent 
crystal field, which controls the occupation of the impu- 
rity orbitals. Fe has ~ 6 d electrons, five of which 
occupy the spin-up d orbitals. The remaining spin-down 
d electron determines the magnetic easy axis (see supple- 
ment for details). If it occupies an orbital with \l z \ > 0, 
one obtains an out-of-plane easy axis, whereas occupa- 
tion of the d z 2 orbital (l z = 0) leads to an in-plane 
easy axis. The resonance in the hybridization function 
of the in-plane <i-orbitals at —2 eV indicates a resonance 
of the Bi2Se3 host at this energy which strongly couples 
to the Fe (d xyi d x 2_ y 2) orbitals. As a consequence, these 
Fe orbitals are pushed upwards within the energy region 
-2eV < E < OeV (c.f. Fig.^c) and (d)) due to level re- 
pulsion between the resonance of the Bi2Se3 host and the 
Fe (d xyi d x 2_ y 2) orbitals. Therefore, the Fe (d xy ,d x 2_ y 2) 
orbitals become depopulated, which facilitates in-plane 
orbital momenta and an in-plane magnetic easy axis (see 
supplement for a detailed analysis in terms of an Ander- 
son impurity model). 

The crystal fields obtained from fitting the spectra of 
CF multiplet calculations to the XAS/XMCD experi- 
ments were D a = 0.19 eV and D r = 0.05 eV. These val- 
ues agree qualitatively with the dynamical crystal fields 
concerning the emptied orbitals as seen by the red curves 
at energies between E = — 2 eV and —1 eV in Fig. [3jc) 
and (d). The fitted spectra are thus in line with the 
picture of dynamic hybridization effects leading to a de- 



population of the Fe (d xy , d x 2_ y 2~) orbitals and a resulting 
in-plane magnetic easy axis. 

In summary, utilizing a combination of experiment 
and theory, we demonstrate that the favorable orienta- 
tion of the magnetic moment of Fe impurities absorbed 
on the Bi2Se3 surface is in-plane with a considerable 
anisotropy but no hysteretic behaviour, i.e. magnetic 
ordering. This has important consequences for the sta- 
bility of the topological Bi 2 Se3 surface state against per- 
turbations by magnetic Fe impurities, which break time- 
reversal symmetry: On average, the exchange field by 
the Fe impurities acts like an in-plane oriented magnetic 
field B on the spins of the surface state Dirac fermions: 
Hd = Vf(a • p) + qiibE • a [28]. From this Dirac Hamil- 
tonian it becomes immediately clear, that an in-plane 
exchange field B on average merely shifts the DP to 
Pd — —gfisB/vf but it does not open a gap, despite 
of the breaking of time-reversal symmetry. This qual- 
itatively agrees with new photoemission spectra which 
do not reveal a gap at the DP [9] [10] as previously ob- 
served [7j. Furthermore, we show that the favored orien- 
tation of the moment can be well understood from the 
interplay of local atomic physics with dynamic hybridiza- 
tion effects. This illustrates the necessity to account 
for both, local Coulomb interaction and dynamical hy- 
bridization effects, when considering ab initio approaches 
for these systems. 
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